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INTRODUCTION 


This Mathematics Student Book was written to enrich your study of mathematics 
by showing you how to use a computer as a modeling device. The computer 
is particularly helpful in quickly performing the repetitious steps of algorithms, 
thus making mathematical investigations easier and more exciting. You will write 
computer programs which will help you to understand the major concepts involved 
in the study of a particular mathematical area. If you become more involved in 
investigating the laws of mathematics, this book will have achieved its aim. 





To use the Student Book for Number Sets, you will need the following: First, you 
should have one year’s background in algebra. Second, the book assumes that you 
already know how to write a program in the BASIC programming language, and 
that you understand programming techniques for inputting data, performing alge- 
braic operations, designating variables, assigning values to variables, looping, and 
printing results. If you do not have this background, you will want to study BASIC 
before attempting this material. Consult the BASIC Manual for the computer you 
use. Last, in order to complete the exercises in this Student Book, you will need 
to have adequate access to a computer (about two hours per week for a terminal 
system). If more time is available, you may be able to experiment further on your 
own, either to improve your programs or to investigate other areas of mathematics 
that interest you. 


Fach section of this book is organized in the same way. First, the mathematical 
concepts needed to complete the exercises are reviewed. References are listed at 
the end of each section in case you want to study these concepts in greater detail. 
Next, each exercise is presented. Finally, an approach is suggested in the Problem 
Analysis and a flow chart is included to illustrate this approach. The suggested 
approach was chosen because it brings out the concepts which are being stressed, 
but the program can sometimes be written more efficiently. Once you have com- 
pleted the exercise by following the logic in the flow chart, you are encouraged to 
rewrite the program using more sophisticated programming techniques. You might 
also want to impose more conditions on the problem to make it more interesting 
to solve. 


There are many different ways to solve one problem by programming. Experi- 
ment and learn as you go. You'll find you are learning something new each time, 
both about your subject matter and about using the computer to solve problems 
and model mathematical algorithms. 
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LIST OF SYMBOLS 


FLOW CHART SYMBOLS 


aren Seis else ai Start or Stop 


Jueesoee ees Defines a process 


<> aa---------- Represents a decision point 


qa iiare acres oe Represents computer input 


arb. 2. oan e Represents computer output 


Used to connect one part of a flow 
chart continued at some other place 
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ALGEBRAIC NOTATION AND EQUIVALENT BASIC LANGUAGE SYMBOLS 





Algebraic Notation BASIC Notation Meaning 
+ + Addition 
- — Subtraction 
- or X ‘ Multiplication 
>or/ / Division 
VJ X SQR(X) Square root of X 
| Y | ABS(Y) Absolute value of Y 
| X | INTCX) Greatest integer less 
than or equal to X 
= = Equals 
# # or <> Does not equal 
— < Less than 
> > Greater than 
<q <= Less than or equal to 
z= >= Greater than or 
equal to 
> = Replaced by 
()orl[ ] ()or[ |] Inclusive brackets or 
parentheses 
A, A(T) Subscripted variable 
A. jj A(I,J) Double subscripted 
variable 
(None) RND(X) Assign a random 


number to the 
variable X 


NUMBER SET DESIGNATIONS 


N — Natural number set Q — Rational number set 
W — Whole number set Z — Irrational number set 
I] — Set of integers R — Real number set 


NUMBERS AND MATHEMATICS 


It may seem obvious to you that numbers are a necessary element of mathematical 
study. As a matter of fact, they have contributed to the development of mathe- 
matics in two important ways: 





First, numbers are used to express quantitative relationships between objects or 
concepts to which number values have been assigned. In this role, they have pro- 
vided the very basis for organized mathematical study. Numbers have served this 
utilitarian role in business, economics, science, engineering, etc. In recent years, 
the computer has made it possible to apply numerical methods to other fields 
such as medicine, science, and social systems. 


Second, number theory (the study of numbers) is an important area of 
mathematics. From the beginning of the use of numbers, man has been curious 
about them. This curiosity first emerged as a form of mysticism ascribed to num- 
bers. But as time went on, the interest in the properties of numbers developed 
into a serious study of their underlying principles. The study of number theory 
was not undertaken for practical reasons, but it has been found to have many 
practical applications. 
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INTEGERS 


“God created the natural numbers, everything else is man’s handiwork.”’ 
—— Leopold Kroeneckes (1823-1891) 


The properties of the set of integers have fascinated mathematicians for hundreds 
of years. The Pythagoreans, a group of mathematicians led by Pythagoras (572- 
501 B.C.), probably did the first serious study of the integers. Their investigation 
included the study of certain series of integers, factors, prime numbers, etc. Euclid 
(300 B.C.) extended the study to include number pairs. 


The computer has made it possible for mathematicians to greatly extend the scope 
of all previous studies of integers. For example, the Pythagoreans investigated 
perfect numbers. Euclid went on to develop a formula, 29-! (2" - 1), that yields a 
perfect number when 2" — 1 is a prime number. Yet by 1876, only twelve perfect 
numbers had been found, mainly because of computational limitations. The 
twelfth perfect number is 2!%6 (2!4” — 1) = 137,438,691,328. You can see why 
further calculations were abandoned until the computer became available. There 
are now 23 known perfect numbers. The 23rd is 2!!! (2!12!3 — 1), which has 
6,751 digits! 


Integers are basic not only to the real number system but to the actual design of 
the digital computer. Therefore, it is appropriate to begin our study of number 
sets with them. 


GREATEST COMMON DIVISOR 


There are many arithmetic operations on the set of integers for which we need to 
find the greatest common divisor. Let’s define the greatest common divisor (g.c.d.) 
in the following way: if dis a common divisor of a set of integers, and every other 
common divisor of that set is a divisor of d, then d is the greatest common divisor 
(g.c.d.). One method for finding the g.c.d. of a set is educated guesswork. 
Start with the smallest integer in the set. If it is a common divisor for all other 
members of the set, it is the g.c.d. Otherwise, factor it and test each of its factors 
until you find one which is a common divisor of the set. If there are no factors 
left greater than the common divisor found, it is your g.c.d. If there are larger 
untested factors, test them one by one until you are left with the g.c.d. If there is 
no common divisor for the set, the integers in the set are said to be relatively prime. 


EXERCISE | — Finding the Greatest Common Divisor 


Write a computer program that will find the greatest common divisor (g.c.d.Jofa 
given set of positive integers. Test your program on the sets G ={91 26,169,286} 
and F = {71,251,149,353}. If the elements of the set are relatively prime, have the 
computer so indicate. 


Integers 
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Problem Analysis 


The flow chart given here uses uneducated guesswork as a method for solving this 
problem. It tests all integers from 1 to some element of the set as common divisors. 
Once you have written the program described by it, make a flow chart of the more 
efficient method described above and alter your program accordingly. 


Macro Flow Chart 


Exercise |. 






INPUT: GIVEN 
SET OF 
INTEGERS 


DETERMINE FIRST 
POSSIBLE 
COMMON DIVISOR 

(c.d.) 












c.d. 
DIVIDE SET OF 
INTEGERS 





NEXT 
POSSIBLE 
c.d. 








SET 
RELATIVELY 
PRIME 


OUTPUT 
g.c.d. 




















“SET RELATIVELY 


EUCLIDEAN ALGORITHM 


Euclid devised an algorithm for finding the greatest common divisor of two posi- 
tive integers. This Euclidean Algorithm, as it is called, is a sequence of steps apply- 
ing the division algorithm. 


The division algorithm tells us that for any integers a and b, there exists integers 
q and r such that a = b(q) + r or in another form, a — b(q) = r. For this dis- 
cussion, we will accept this division algorithm without proof. Given two integers, 
the Euclidean Algorithm substitutes the integers for a and b in the division al- 
gorithm. The variable q is given the value of the greatest integer for which b(q) <a 
and then the value of r is found. ais then replaced with b and b withr. This pro- 
cedure is repeated until r equals zero. 


The final value of b is the g.c.d. of the original values of a and b. To continue to 
check the elements of the set, we use the g.c.d. just found for a and the next 
element in the set for b. We continue the procedure until all elements in the set 
have been considered. Now let’s apply the algorithm to the set below: 


A = {81,54,1458,36 } 


To find the g.c.d. of 81 and 54, the first two elements of A, we follow our 
procedure: 


a<81,b< 54 gives us 81 — 54(1) = 27. 


The next step gives us a < 54, b « 27 and 54 — 27(2) = 0. Therefore, the g.c.d. is 
this final value of b, or 27. 


The next element in the set is 1458. We start our calculation with a equal to the 
g.c.d. already found, 27 in this case, and b equal to 1458: 


27 — 1458(0) = 27 

a<- 1458, b< 27 

1458 — 27(54) =0 
Therefore: The g.c.d. of 27 and 1458 is 27. 


a «+ 27, b < 36 (last integer in set) 
27 — 36(0) = 27 

a<« 36,b< 27 

36 — 27(1) =9 

a<-27,b<9 

27 —9(3) =0 


Therefore: The g.c.d. of the integers in Set A is 9. 
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EXERCISE 2 — Using the Euclidean: Algorithm 


Write acomputer program which uses the Euclidean Algorithm to find the greatest 
common divisor ofa set of integers. Test your program on Sets G and F of Exercise 1. 


Problem Analysis 


Once you understand the above example, you should have no trouble completing 
this exercise. 


Macro Flow Chart 


Exercise 2. 


START 






INPUT: GIVEN 
SET OF 
INTEGERS 


APPLY ALGORITHM 
TO FIRST PAIR OF 
INTEGERS, a AND b, TO 
DETERMINE THE g.c.d. 














ALL 
INTEGERS 
BEEN 
CONSIDERED 


=eS OUTPUT: 


g.c.d. 





NO 


APPLY ALGORITHM STOP 
TO g.c.d. AND NEXT 


INTEGER OF SET TO 
DETERMINE NEW g.c.d. 





10 


PRIME NUMBERS 


A prime number is a positive integer which has no divisor other than itself and the 
number one. The most obvious method for determining primeness is simply trial- 
and-error. We’ll use this method for the next exercise. 


EXERCISE 3 — Prime Numbers 


(a) Write a computer program that will determine all prime numbers between I 
and 500. 


(b) After you have a working program for Part (a), modify your program so it 
will count the number of prime numbers in the intervals 1-100, 1-200, . . . 
1-500. 


(c) Adjust the program for Part (a) to count and print twin primes from I to n 
asn < 100, 200,... 500. Twin primes are consecutive primes. For example, 
3 and 5 are twin primes. 


(d) Adjust program 3(a) to have the computer print out ratios n/m and n/logen 
forn< 107, 10°, 10%, 10°, 10°, where m is the number of primes less than n. 


Problem Analysis 


The flow chart clearly illustrates the solution for Part (a). No flow charts are in- 
cluded for Parts (b), (c), and (d) since they are merely modifications of Part (a). 
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Macro Flow Chart 


Exercise 3(a) 


n <—FIRST NUMBER 


TO BE TESTED FOR 
PRIMENESS 


















d < FIRST POSSIBLE 
DIVISOR OF n TO 
BE CONSIDERED 





n<— NEXT NUMBER 
TO BE TESTED 











ALL 







d MBER 
A DIVISOR pers > 
OF n CONSIDERED 








ALL STOP 











NECESSARY 
d < NEXT POSSIBLE DIVISORS 
DIVISOR OF n OF n BEEN 
CONSIDERED 





OUTPUT: 
VALUE OF n 


2 


EXERCISE 4 — Determining Prime Numbers by “Sieve of Eratosthenes’’ Method 


Write a computer program to determine all prime numbers from 1 to some number 
n using the “Sieve of Eratosthenes’”’ approach. 


Problem Analysis 


The sieve approach involves picking the first prime number in the set and then 
eliminating all multiples of that prime. This is repeated until only prime numbers 
remain. For example, we’ll use this approach to find all prime numbers between | 
and 36. We’ve shown the set under consideration and the output generated after 
each step. 


Input: 2 3 4 5 6 7 8 
9 10 11 12 13 14 15 

16 17 18 19 20 21 22 

23 24 25 26 2] 28 29 

30 31 32 33 34 35 36 


First Output: 2 3 5 7 

(Eliminating 9 11 13 i) 

multiples of 17 19 21 

2 except 2) 23 25 DA 29 
31 33 35 

Second Output: 2 3 5 7 

(Eliminating 11 13 

multiples of 17 19 

3 except 3) 23 25 29 
3] 35 

Third Output: 2 3 5 7 

(Eliminating 11 13 

multiples of 17 19 

5 except 5) 23 29 
3] 


The steps illustrated lead us to construct the following flow chart. 
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Macro Flow Chart 


Exercise 4. 










INPUT 
GIVEN SET OF 
INTEGERS FROM 
1TOn 










OUTPUT 
GIVEN SET IN 
MATRIX FORM 

(OMIT THE NUMBER 

ONE FROM THE SET) 








ELIMINATE 
MULTIPLES 
OF p 


OUTPUT 
REMAINING 


INTEGERS 
IN THE SET 


NEXT 
p< PRIME 
NUMBER 





YES 
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PERFECT NUMBERS 


A number, n, is called a perfect number when the sum of all its divisors (including 
1 but excluding itself) is equal to n. 


EXERCISE 5 — Determining Perfect Numbers 


Write a program to find all perfect numbers less than or equal to 500. Use the fol- 
lowing output format for your results: 


Sito oo etait 


where n is the perfect number and f,,f,,... fare the factors ofn such that f,< n. 
Problem Analysis 


This problem is solved by merely testing each number from 1 to 500 to see which 
numbers meet the condition defined above. 


Macro Flow Chart 


Exercise 5. 


SELECT FIRST 
NUMBER, n 
TO BE TESTED 


DETERMINE ALL 
nwa nt FACTORS OF n 
LESS THAN n 


SUM 
OF FACTORS 
EQUALS 
n 








YES OUTPUT 


FACTORS, n 





NO 
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CAN YOU FIND THE NEXT ONE? 





According to Tecnica Education Corporation Newsletter, December 
1970, eighteen-year-old Roy Ferguson of Dallas, Texas recently devel- 


oped a computer program that determined the 21st perfect number. 
The number has 5834 digits. 


PYTHAGOREAN TRIPLES 


Positive integer solutions to the Pythagorean equation z” = x2 + y2 can be found 
by using the relationships (1) x = m? — n?, (2) y = 2mn, and (3) z = m? + n?2 
where m and n are positive integers with m > n. These solutions are called Py thag- 
orean triples. 


EXERCISE 6 — Finding Pythagorean Triples 


(a) Using the relationships defined above write a computer program to find 100 
Pythagorean triples satisfying z? =x* +y?, with x, y, and z positive integers. 


(b) Adjust your program so that it prints only primitive triples, ie., triples that 
have no common divisors greater than 1. 


Problem Analysis 


(a) Our task is to find 100 ordered pairs, (m,n,), with m >n, so that the values 
of x, y, and z obtained by using the given relationships satisfy the condition 
z* i x? — y*. 


Since m >n > 0 we begin with 


First (m,n) = (2,1) 
Next (m,n) = (3,1) 
Next (m,n) = (3,2) 
Next (m,n) = (4,1) 


The flow chart is clear from this point on. 


(b) Your program will only need to determine that x and y are relatively prime. 
Can you verify that this will assure that x, y, and z have no common divisors 
except 1? 
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Exercise 6(a) 


START 


SELECT FIRST 
POSSIBLE ORDERED 
PAIR (m,n) 





YES 












HAVE 
NEXT PAIR NO 100 TRIPLES 
(m,n) BEEN 


FOUND 


YES 


STOP 
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RATIONAL NUMBERS 


A rational number is a number which can be expressed as a/b, where a and b are 
integers with b # 0. The development of rational numbers was prompted by the 
need to measure quantities which were smaller than whole units. 


Early man was slow to develop a number system that included numbers represent- 
ing “aliquot parts” (fractions) of a unit. Fora long time, it was considered adequate 
to measure quantities to the nearest unit. The Egyptians treated only fractions 
with a numerator of 1, and the Babylonians did not master fractional numbers 
until 2000 B.C. 


Abstractly, the need for rational numbers is justified by the need to be able to 
solve equations of the form ax + b = c, where a, b, and c are integers anda #0. 


We will use the computer to investigate the property of denumerability of rational 
numbers. The set of rational numbers is said to be denumerable because it is 
“countable,” i.e., it can be shown to be in one-to-one correspondence with the 
natural numbers. Figure 1 illustrates one way the set of positive rational numbers 
can be denumerated. This method is thoroughly described in Courant’s, What 


Is Mathematics. 


1 /1——2/] 3/1 a 
ae a ae - 
Fc 2, a s 
ou “ 5/4 6/4 
a 3/5 4/5 5/5 6/5 
1/6 2/6 3/6 4/6 5/6 6/6 


Figure 1. Denumerating the Set of Positive Rational Numbers 
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If you follow the path of the arrow in Figure 1, it is obvious that the positive 
rational numbers can be denumerated. If all numbers on the path are denumerated, 
some numbers will be listed more than once (2/1 and 4/2, for example). 


EXERCISE 7 — Denumeration of the Set of Positive Rational Numbers 


(a) Write a computer program that will print the first 100 rational numbers of 
the list generated as described above. 


(b) Alter your program so the output includes zero and the negative rational 
numbers. 


Problem Analysis 


Our problem, then, is to develop an algorithm that will generate the above sequence 
of numbers, printing and counting the rational numbers 1/1, 1/2, 2/1, 1/3, ... 
until 100 distinct rational numbers have been printed. 


Let’s examine the numbers along each diagonal of the path: 


(a) 1/1 

(b) 2/1,1/2 

(c) 1/3, 2/2, 3/1 

(d) 4/1, 3/2, 2/3, 1/4 

(e) 1/5, 2/4, 3/3, 4/2, 5/1 

(f) 6/1, 5/2, 4/3, 3/4, 2/5, 1/6 


and so on. 
Note that the numbers increase or decrease within each segment (diagonal) of the 


sequence. Changing the order of the numbers within a segment would not void the 
process but would simplify the programming of the process. 


You should realize that the output of this algorithm will not list the rational num- 
bers in order of magnitude, but it will illustrate the ability to produce the whole 
set of rational numbers, if the computer runs indefinitely. 


20 





NUMBER SETS 





Macro Flow Chart 


Exercise 7(a). 







START 


GENERATE FIRST 
RATIONAL NUMBER, on, 
OF THE SEQUENCE 












n 
EQUAL TOA 
PREVIOUSLY 
STORED n 


100 
RATIONAL 
NUMBERS 
PRINTED 












GENERATE NEXT 
RATIONAL 
NUMBER n 







EXERCISE 8 — Forming Rational Numbers, Using All Integers in a Given Real 
Number Interval 


(a) Given a set of integers, W= {x |-75xS 7}, have the computer print and 
count all meaningful rational expressions that can be formed from set W. 


(b) Adjust your program so the computer will print only those expressions that 
do not duplicate values of expressions already printed. Also, have the com- 
puter output the decimal approximations of each number printed. 


(c) Adjust the program from part (b) to print all expressions in reduced form. 
Problem Analysis 
(a) If we were doing this assignment by hand, one approach might be to write: 


~7/-7, -7/-6, -7/-5 .. . -7/-1, -6/-7, -6/-6, -6/-5, ... 
0/-7, 0/-6, 0/-5,...0/7... 7/7. 


The flow chart below explains the principle behind this approach. 


Rational Numbers 


21 





MATHEMATICS 


Hewlett-Packard Computer Curriculum 


Macro Flow Chart 


Exercise 8(a). 












INPUT 
GIVEN SET 
OF INTEGERS 


a < FIRST 
POSSIBLE 
NUME RATOR 


b «+ FIRST 
POSSIBLE 
DENOMINATOR 











ae 

*) 

/ \ 
a/b YES 

MEANINGFUL oor 

EXPRESSION ¢,a''/""b 














ALL 
POSSIBLE 
DENOMINATORS 
USED 


NO 
b <— NEXT 


DENOMINATOR 


YES 














ALL 
POSSIBLE 

NUMERATORS 
USED 


NO 
a<«< NEXT 


NUMERATOR 


YES 


STOP 


NOTES: 


1. The section in the circle marked (*) will be replaced by another section to 
meet the conditions of Part (b). 

2. a/b indicates quotient of a divided by b. 

3.  a‘‘/’b represents the rational expression of a “‘over”’ b. 


Le 








Part (b). 


Substitute the following section for the circle in the flow chart for Part (a). 











a/b 
EQUIVALENT 
TO ANY 
PREVIOUSLY PRINTED 
EXPRESSION 






OUTPUT 
c,a‘’/""b,a/b 


A flow chart section will be inserted at (*) to satisfy the conditions of Part (c). 


Part (c) 


Insert the following section as indicated into flow chart of Part (b). 










REDUCE a’’/"b 
TO LOWEST TERMS 


EXERCISE 9 — Rational Numbers Formed from a Random Set of Integers 


Given a set of random integers, A = {—] 6, 3,45, 0, 24, -6, 7,19, 38}, have the 
computer form all possible unique rational numbers, separating them into three 
categories: 


(a) -1 <plq <0 


(b) O<plq <1 
(c) | p/q|>1,pandqe A. 


Have the computer output results in column form, printing the decimal approxt- 
mation of each rational number. 


Problem Analysis 


This is a variation of the general case considered in Exercise 8. The procedure is 
the same. 
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Macro Flow Chart 


Exercise 9. 
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START 


PRINT TABLE 
HEADINGS 


p< FIRST ELEMENT OFA 
q~ FIRST ELEMENT OFA 






p/q 
RATIONAL 
NUMBER 











OUTPUT 

























DETERMINE p’/"q 
TO EXPRESSION PROPER AND DECIMAL 
ALREADY CATEGORY APPROX. IN 
PRINTED FOR p/q PROPER COL. 









ALL 
POSSIBLE 
VALUES FOR 
q USED 


p< NEXT ELEMENT OFA 
q < FIRST ELEMENT OFA 
ALL 


POSSIBLE 
VALUES FOR 
p USED 








q < NEXT 
ELEMENT 
OFA 












YES 


STOP 
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IRRATIONAL NUMBERS 





Early mathematicians were convinced that all points on the number line are used 
up when a point is assigned to each rational number. About 500 B.C., the Pythag- 
oreans proved that there are points on the line which do not correspond to rational 
numbers. In particular, they proved that no rational number corresponds to the 
measure of the diagonal of a unit square. A new category of numbers was invented, 
called irrational numbers, to represent known values which could not be expressed 
as rational numbers. Irrational numbers find widespread application today. 


NEWTON’S METHOD FOR FINDING SQUARE ROOT 


One method of determining a decimal approximation to the square root of a num- 
beris known as Newton’s method. It is also called the successive approximation or 
averaging method. Whatever we call it, we start by making a guess at the square 
root of a number, n. We'll call this guess g,. If g, >4/n, you can see that n/g, < 
Jn. Or, if g, </n, then n/g, >/n. We'll assume that g, > /n. A logical sec- 
ond guess, g,, would be the average of g, and n/g,. We continue the process by 
picking £4 equal to the average of g, and n/g, and so on until our guess, g,, is 
within a given tolerance of the actual value of the square root of n. Figure 2 shows 
the first two steps on the number line. 


n/g, Jn 81 


<———<$ $$ $$ rs 
g, = 1/2(g, + n/g,) 


nig, fn & 


eo mn 0— 
g, = 1/2(g, + n/g,) 


Figure 2. A Number Line Representation of Newton’s Method 


EXERCISE 10 — Rational Approximation of \/ n, n rational > 0 


(a) Using Newton’s method, write a computer program that will find a decimal 
approximation to the square root of any positive rational number accurate to 
within .001 of the actual value. Test your program on several numbers. Com- 
pare the results with a square root table. 


(b) Adjust the program of Part (a) to find the cube root of a positive rational 
number. Will your program find cube roots of negative numbers? 


(c) Now, you might wonder, since we are taking a cube root, wouldn ta good 
guess be 1/3(g;_, + n/g;-;7)? Try it and see what happens. 


Irrational Numbers 
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(d) We will not be able to show the development here, but it can be shown that 
g; — 1/3(2g;., + n/g;.;*) is another good substitution to use to find the cube 
root of n. Use it in your program. Compare your results with those obtained 
in the first cube root program. Do some research on why the algorithm will 
produce cube roots. 


Problem Analysis 
(a) The flow chart clearly follows the procedure discussed above. 


(b) This part of the exercise requires some modification to our procedure. We 
start with a guess, g,, at the actual value of on. Let’s assume g) ></ Ui, 
which means that n/g,? <<) ie Of course, if our guess is accurate, g, r= a 
and n/g, — g,. But guesses aren't very often that accurate. 


We zero in on the actual value by making g, equal to the average of g, and 
n/g,7, and so on. In general, then, g, = 1/2(g_, + n/g? | ). Substitute this 
formula into your program for Part (a). 


Macro Flow Chart 


Exercise 10(a). 






INPUT 
ANY POSITIVE 
RATIONAL 

NUMBER, n 














INPUT 
g << INITIAL 
GUESS AT \/n 









IS g 
WITHIN THE .001 
TOLERANCE 


YES OUTPUT 
g 










g< THE 
AVERAGE OF 
g AND n/g 
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INCREMENTING METHOD 


A less sophisticated method for finding the nth root of a number is the increment- 
ing method. This method demonstrates the power of the computer to perform a 
systematic search for a desired value within a certain tolerance. 


Suppose we wish to find 3/76 within .001 of the actual value. We make a guess g, 
an integral value. Assume we guess g = 2, then we test g* and find 2? < 76. Then 
we increment g so g= 3. g? is still less than 76 so we let g = 4, then 5 until we find 
5? > 76. Therefore, we know that 3/76 is between 4 and 5. We set g back to 4 
and increment g by .1 until g? > 76. Again, we reset g to the previous value so that 
g> < 76, and repeat the process with the next increment. We continue this pro- 
cedure until we reach the desired degree of precision. We find that 4.235 <</76 
< 4.236, and either 4.235 or 4.236 is an acceptable solution. 


We might find that the initial guess, g, was greater than 3/ 76. In that case, the same 
procedure is used after the value of g is decremented until g<./ 76. 


EXERCISE 11 — Incrementing Method for Finding \/ n, n rational 


Write a computer program that will find a number g <rth root of a positive 
rational number n, so that n—g' <t where t is a given tolerance. 


The output of the problem will be much more interesting if the program is 
written so that each tested value of g is printed. In that way, you can examine 
the searching process used by the computer prior to the final output. 


Apply your program to the following cases: 
1. cube root of 16, t=.001 


2. fourth root of 81, t = .001 
3. square root of .0234, t= .00001 


Problem Analysis 


This problem is a very straightforward application of the incrementing method. 
When incrementing make the new increment equal to one tenth of the previous 
increment. 
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Macro Flow Chart 


Exercise 11. 
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START 













INPUT 
n, NUMBER WHOSE ROOT IS TO BE DETERMINED 
r, ROOT DESIRED 

g, INITIAL GUESS FOR r*4 ROOT OF n 

t, TOLERANCE TO WHICH ROOT IS DESIRED 


K<-1 




















DECREASE GUESS g 
BY INCREMENTS OF k 


UNTIL gi < n 


K< K/10 


INCREASE g 
BY ONE 
INCREMENT, k 


g 
WITHIN 
THE REQUIRED 
TOLERANCE 


STOP 
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Il, THE MOST FAMOUS [IRRATIONAL NUMBER 


The numberIlis undoubtedly the most famous irrational number. It has fascinated 
man since he first discovered that it existed. Since then, mathematicians have made 
it the subject of innumerable studies, most of them trying to determine one of the 
following: 1) the exact value of II, and 2) the distribution of its digits. In fact, for 
centuries an unbelievable amount of effort was expended in computing I] to as 
many decimal places as possible. 


The value of II was first established in ancient Oriental civilizations as the ratio of 
the diameter to the circumference of a circle. The earliest value given to it was 3, 
according to Biblical reports (1 Kings 7:23; I] Chron. 4:2). About 240 B.C., Archi- 
medes computed it to be between 223/71! and 22/7 or approximately equal to 
3.14. By 150 A.D., it had been calculated to be 3.1416. 


During the next 1500 years, many mathematicians struggled with computing II to 
more and more decimal places. All used some approach related to the circle. The 
last major effort using this approach was made in 1630 by Grienberger, who was 
able to compute II to 39 decimal places. 


By the end of the 18th century, mathematicians had demonstrated that Ilis not 
the exclusive property of circles. Various series were developed which can be used 
to compute II (reference Eves’ and Davis’ books for details). Using such formulas, 
the value of II was computed to more and more places until, in 1873, William 
Shanks of England computed it to 707 places. It took him 15 years to complete 
his calculations. 


In 1940, Edward Kasner and James Newman, in their book Mathematics and the 
Imagination, stated: ‘even today it would require 10 years of calculation to deter- 
mine II to 1,000 places.” Yet in 1949, by use of the electronic calculator ENIAC, 
the Army Ballistic Laboratories in Aberdeen, Maryland produced II to 2035 places 
in 70 hours. Today, II has been computed to more than 100,000 decimal places. 


Two interesting approaches to computing II are found in the area of probability. 
The first is the Buffon needle problem which we will not attempt to explain here. 
Refer to Gamow’s book for an explanation of this problem. The other approach is 
based on a relationship of randomly chosen pairs of positive integers. We will 
investigate this second approach. 


The relationship is as follows: the probability that two randomly chosen positive 
integers are relatively prime is 6/II*. This means that if we consider an infinite 
number of pairs of randomly chosen numbers and tabulate those pairs that are 
relatively prime as successes, s, and those that are not as failures, f, then s/(s + f) = 
6/Tl?. Solving this equation for N, we get IL=/ (6s + 6f)/s. 
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EXERCISE 12 — Computing I by a Probability Method 


Write a computer program that will consider n pairs of random positive integers 
and, after each pair has been tested for relative primeness, evaluate and print 
J (6s + 6f)/s. Run your program for n = 100, 200 and 500. 


Problem Analysis 


You will be able to compute 7 to a surprising degree of precision by using only 
random numbers selected from the set of integers | to 1000. Of course, you can 
use numbers larger than 1000 if you desire. 


Macro Flow Chart 


Exercise 12. 
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NEXT RANDOM 
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